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Abstract. Ultra-high-energy (E>100 Tev) Extensive Air Showers (EASs) have 
been monitored for a period of five years (1997 - 2003), using a small array of 
scintillator detectors in Tehran, Iran. The data have been analyzed to take in 
to account of the dependence of source counts with respect to the zenith angle. 
During a calendar year different sources come in the field of view of the detector at 
varying zenith angles. Because of varying thickness of the overlaying atmosphere, 
the shower count rate is extremely dependent on zenith angle which have been 
carefully analyzed over time. fBahmanaba di et al. 20 02). High energy gamma-ray 
sources from EGRET third catalogue where observed and the data were analyzed 
using an excess method. A number of upper limits for a number of EGRET sources 
were obtained, including 6 AGNs or probably AGNs and 4 unidentified sources. 

Key words. EGRET sources, Extensive Air Showers (EASs), Gamma- Ray sources 



1. Introduction 

EGRET instrument on-board Compton Gamma Ray Observatory (CGRO) has de- 
tected both diffuse and discrete gamma-ray emission. The diffuse emission is both 
galactic (Hun ter et al. 1997j) and extra-galactic in nature IjSreekumar et al. 1998(1 . 
EGRET has detected about 271 high energy (>100 MeV) gamma-ray sources 
IjHartman et al. 1999(1 . Besides AGNs these sources include 170 sources that are not 
identified conclusively with unique counterparts in other wavelengths. Two third of 
these EGRET unidentified (EUI) sources lie close to the galactic plane-potential 
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counterparts ( |Bhattacharya et al. 2003| ) for these include young pulsars, young radio 
quiet pulsars IjTorres et al. 20011 ID'Amico et al. 20011 |Zhang et al. 2000| ), Wolf Rayet 
(WR) stars, Of stars, OB associations l|Romero et al. 1999)) . Super Nova Remnants 
(SNRs) l)Oombi et al. 20011 |Case fc Bhattacharya 19981 ISturner fc Dernier 19951 and 
other types of sources. 

Some other faint sources are in the mid-latitude region suggested to be associ- 
ated with the Gould Belt (Gch rels et al. 2000J) . which underwent an intense star for- 
mation period about sixty million years ago IjGrainer 20001 |Harding fc Zhang 2001| ). 
High latitude sources which are about 50, might be galactic gamma-ray halo sources 
l)Dixon et al. 1998)) or unidentified sources are thought to be extra-galactic. These extra- 
galactic EUI sources contain Blazars and Active Galactic Nuclei (AGNs), galaxy clusters 
< |Colafrancesco, S., 2002| ), BL Lacerta objects IjTorres et al. 200*3)1 and other types. 
Whether the EGRET sources have emission in higher energies, is an interesting ques- 
tion l|Lamb fc Macomb 1997)1 . Gamma-ray with energies about 100 TeV and more enter- 
ing the earth atmosphere, produce Extensive Air Showers (EASs) I jGaisser, T.K., 1990| ) 
which could be observed by the detection of the secondary particles of the showers on the 
ground level IjBahmanabadi et al. 19 98) . Previous attempts have been reported by other 
EAS arrays l)Amenomori et al. 20021 12TTO1 IBorione et al. 19971 lAlexandreas et al. 19931 
|McKay et al. 1993| ). 

This paper reports the results of a small particle detector array located at the 
Sharif University of Technology in Tehran. This small array is a prototype for 
a larger EAS array to be built at an altitude of 2600 m (= 756 g/cm 2 ) at 
ALBORZ Observatory (AstrophysicaL oBservatory for cOsmic Radiation on alborZ) 



(see http://sina.sharif.edu/~observatory/! near Tehran. The prototype installed on the 
roof of physics department of Sharif University of Technology in Tehran, 1200 m (= 890 
g/cm 2 ), 35.72°iV and 51. 33°E. In this work we present the observational results of 10 
EGRET third catalogue sources, we describe the experimental setup in Section 2. the 
data analysis in Section 3, the results in Section 4. Section 5 is devoted to a discussion 
of the results. 



2. Experimental arrangements 

The array is constructed of 4 slab plastic scintillators (100 x 100 x 2 cm 3 ) as a square in 
Tehran (35° 43' , 51° 20'), Iran, with The elevation 1200 m over sea level (890 g/cm 2 ) 
which is shown in Fig.^ All of the scintillators are on a flat level surface. Each scintillator 
is housed in a pyramidal steel box with height of 15 cm. The interior surface of each box 
is coated with white paint, (Bahmanab adi et al. 1998)) and a 5 cm diameter PMT(EMI 
9813KB) is placed at the vertex of the pyramidal box. Fig.^shows a schematic diagram 
of the array and its electronic circuit to log each EAS event. After passing of at least one 
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particle from a detector the PMT creates a signal with a pulse height which is related to 
the direction, number of the passed particles, and location of the crossed particles in the 
scintillator. The output signals from the PMTs are amplified in a one stage amplification 
(xlO) with an 8-fold fast amplifier (CAEN N412), and then transfer to an 8-fold fast 
discriminator (CAEN N413A) which is operated in a fixed level of 20mV one by one. 
The threshold of each discriminator is set at the separation point between the signal and 
background noise levels. Each discriminator has two outputs, one of them is connected to 
a coincidence logic unit (CAEN N455) as trigger condition. Trigger condition is satisfied 
when at least one charged particle passes through each of the four detectors within a time 
window of 150ns. The other discriminator output is connected to a Time to Amplitude 
Converter (TAC)(EG&G ORTEC 566) which are set to a full scale of 200ns (maximum 
time difference between each two scintillators which is acceptable) . The outputs of the 
No. 4 scintillator was connected to start input of TAC1 whereas the output of No. 2 was 
connected to start inputs of TAC2 and TAC3. The Output of the scintillator No. 3 was 
connected to the stop input of TAC2 and No. 1 was connected to stop inputs of both TAC1 
and TAC3. Then the outputs of these three TACs were fed into a multi parameter Multi 
Channel Analyzer (MCA)(KIAN AFROUZ Inc.) via an Analogue to Digital Converter 
(ADC)(KIAN AFROUZ Inc.) unit. 

When all of the scintillators have coincidence pulses, these TACs are trigged by logic unit 
and 3 time lags between the output signals of PMTs (4,1), (2,3) and (2,1) are read out 
by a computer as parameters 1 to 3. So with this procedure an EAS event is logged. 
Two different experimental configurations were used by the experimental set up. All of 
the experimental set up were identical in the first (El) and the second (E2) experimental 
configurations except for the size of the array. In El the size was 8.75 m x 8.75 m and 
in E2 the size was 11.30 m x 11.30 m. 



3. Data Analysis 

The logged time lags between the scintillators and Greenwich Mean Time (GMT) of 
each EAS event were recorded as raw data. We synchronized our computer to GMT 



(see http://www.timeanddate.coml. Our electronic has record capability of 18.2 times 
per second or equivalently each 0.055 seconds one record will be stored regardless of 
the existence or non existence of EAS events. If an EAS event occurs, its three time 
lags will be recorded and if it does not occur 'zero' will be recorded. Therefore with 
the starting time of each experiment and counting of these zero and non zero records 
we will obtain GMT time of each EAS event. Our detected EAS events are a mixture 
of cosmic-ray events and gamma-ray events. In El total number of EAS events was 
53,907 and duration of the experiment was 501,460 seconds. So the mean event rate of 
the first experiment was 0.1075 events per second. The distribution of the time between 
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successive events has a good agreement with an exponential function, indicating that the 
event sampling is completely random ijBahmana badi et al. 2003^1 . In E2 total number of 
events was 173,765 and duration of the second experiment was 2,902,857 seconds. So its 
mean event rate was 0.05986 events per second. 

We refined the data for separation of acceptable events. Events are acceptable if there 
be a good coincidence between the four scintillator pulses. We omitted the events with 
zenith angles more than 60°. Therefore after the separation we obtained smaller data 
sets of 46,334 and 120,331 for El and E2 respectively. Since we can not determine the 
energy of the showers on an event by event basis, we estimate our lower energy threshold 
by comparing our event rate to a cosmic-ray integral spectrum (jBorione et al. 1997(1 

J(E) = 2.78 x \Q- 5 E- 2 - 22 + 9.66 x lO^iT 1 - 62 - 1.94 x 10~ 12 40 < E < 5000 TeV (1) 

The obtained lower energy limits were 39 Tev in El and 54 TeV in E2. The calculated 
mean energies were 94 and 132 TeV in El and E2 respectively. If the well-known Hillas 
spectrum ( |Gaisser, T.K., 1990| ) 

„, „s m particle E , „ 

f >£ ^2xl0- 10 ^ x — — — " 7 2 

cm/ s sr 1000 J ev 

is used the lower limits will be 40 TeV and 60 Tev. Since the distribution of cosmic- 
ray events within the array in these energy ranges are homogeneous and isotropic, we 
used an excess method l|Amenomori et al. 20021 12000|l to find signature of EGRET third 
catalogue gamma-ray sources. This method was used for both El and E2. 
The complete analysis procedure is itemized as follows : 

— The calculation of local coordinates; zenith and azimuth angles of each EAS event 
(z, ip) were calculated using a least square method by logged time lags and coordinates 
of the scintillators. 

— The local angle distributions of the EAS events were investigated to understand the 
general behaviours of these EAS events. 

— The calculation of equatorial coordinates (RA,Dec) of each EAS event using its local 
coordinates, the GMT of the event and geographical latitude of the array. Then we 
calculated galactic coordinates (l,b) of each EAS event from its equatorial coordinates 
using epoch J2000. 

— The estimation of angular errors in galactic coordinates of investigated EGRET 
sources by error factors of the array. 

— The simulation of a homogeneous distribution of EAS events to investigate cosmic-ray 
EAS events. This simulation incorporated all known parameters of the experiment. 

— The investigation of the statistical significance of random sources and the signifi- 
cance of sources from the third EGRET catalogue using the method of Li & Ma 
IjLi fc Ma 1983|l and find the best location for EGRET sources in the TeV range. 



M. Khakian et al.: Observation of EGRET Gamma-ray sources ... 5 
3.1. Calculation of Local coordinates of each EAS event 

The local coordinates are zenith (z) and azimuth (</?). We used the least square method 
( |Mitsui, K., et al. 1990| ) to calculate z and ip. It is assumed that the shower front could 
be approximated by a plane. So we obtain, 
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D j = Dj — D =.T ji+2/ jj and t j = tj — t are the coordinate vector and the time 
lag of ]th scintillator with respect to the reference one and c is the velocity of light. 
A zenith angle cut off 60° is implemented to enhance significance. 

3.2. Angular distribution of the EAS events 

Fig. |2Ja) shows the azimuthal angle distribution of the EAS events which is nearly 
isotropic. A slight North-South anisotropy is observed which is attributed to the 
geomagnetic field. We fitted this distribution with a harmonic function as follow : 
( Bah manabadi et al. 2 002 ) 

f((j)) = A<f, + cos(0 - ipi) + Ccj, cos(20 - (p 2 ) (6) 

where A^,, B^, and are respectively 14516, 1270 and 184. <p\ and </?2 are phase con- 
stants which are respectively 32° and 55°. 

Since thickness of the atmosphere increases quickly with increasing zenith angle z 
( |Gaisser, T.K., 1990| ), the number of EAS events is strongly related to the z value, as 
shown in Fig. HJb) . 

These distributions were studied separately for the two experimental configuration El 
and E2. The shower rate in E2 is less than El because of the larger size of the ar- 
ray in El. However the zenith angle distributions in El and E2 are very similar. 
The differential zenith angle distributions of these data sets are fitted to the function 
dN = A z SinzCos n zdz with a very good agreement for both El and E2 which from 
0° to 50° A z = 95358 and n = 5.85 and from 50° to 60° A z = 90189 and n = 5.00. 
From another view the mean value of zenith angle (z) is 26.2746 and 26.4625 in El and 
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E2 respectively. Since the results of the two experimental configurations, are in a good 
agreement with one another and the excess is important for us. Therefore we added the 
two data sets to obtain a larger data set with lower energy threshold of El which is 39 
TeV. 



3.3. Calculation of equatorial and galactic coordinates of each EAS event 

The equatorial coordinates (RA,Dec) are obtained from calculated local co- 
ordinates {z,tp), GMT of each EAS event and geographical latitude of the 
array. In this step the transformation relations (see http: / / aanda.u-strasbg.fr| 



|Roy, A.E., fc Clarke, D.| ), and the local sidereal time of the starting point of the ex- 



periment (see http://tycho.usno.navy.mil/sidereal.htmll were used 



Then galactic coordinates (l,b) of each EAS event are obtained from the calcu- 
lated equatorial coordinates, based on the galactic coordinate standard of year 2000 
(see http://aanda.u-strasbg.frl. Fig. shows the distribution of our data in galactic co- 
ordinates. 



3.4. Error estimation of investigated sources in galactic coordinates 

For the coordinates calculations of each EAS event in galactic coordinates we have to 
know estimated errors in these coordinates. These errors are due to experimental error 
factors, which contain uncertainties in time and coordinates of each logged EAS event. 
The defined distance between two scintillators was centre to centre and the size of the 
scintillators were (100 x 100 x 2 cm 3 ). Meanwhile the accuracy of coordinates of each 
scintillator is measured within a few centimeters. So error in measurement of coordinates 
of secondary particles of each EAS event is Ad = 1 m. 

The errors in time measurement of each EAS event are due to the front thickness of 
the secondary particles, electronics errors and error in GMT logging. The error due 
to the first two factors was At = 2ns ( Bahmanabad i"et al. 2002)l . The error in GMT 
logged time of each EAS event was AT = 0.07s which is due to recording rate and the 
synchronizing of the computer. These errors make uncertainties in galactic coordinates 
of the investigated sources by the array. 
The following errors were calculated : 



— The errors in local, equatorial and galactic coordinates of each EAS event. 

— The observational angular error of each source. 

— The mean and standard deviation of these error angles. We calculated these steps for 
more than 1000 random sources which are in the Field Of View (FOV) of the array. 
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This calculation was carried out for El and E2 separately and was weighted by their 
refined EAS events. 



From geometry of Fig. ^ we can drive : 

sinz sin 93 = -(t 2 — t±) (7) 

sinzcos^ = -(t 3 — i 2 ) (8) 
a 

In these relations Us are logged times of an EAS event in i t h component of the array and 
d is side of the square array. 

The errors in zenith and azimuth angles were obtained by differentiating from eqs. (5) 
and (6) : 

Az 2 = A 2 A / cos 2 z + B^tim 2 z (9) 

Atp 2 = A\/ sin 2 z (10) 

where A a = 2cAt/d and Ba = Ad/d. The errors in equatorial and galactic coordinates 
were calculated from differentials of Dec(z,<£>), RA(z,y, T), b(RA,Dec) and l(RA,Dec). 
If y is a generic function of parameters u, v and T, then : 

y = y{u,v,T) (II) 



|A,,| : : \/(^) 2 Au 2 + (^) 2 A V 2 + (^) 2 AT 2 (12) 



Where dy/dT = 2a (a = 1.00273790935) (see |http: / /tycho.usno.navy.mil/sidereal.html I 
for the calculation of RA(z, <p,T) and dy/dT = for calculation of Dec(z,^j), b(RA,Dec) 
and l(RA,Dec). 

The error on the observed solid angle of each source is Afl = cosbAbAl and the equivalent 



error on the angular radius is r e = y/ AQ/ir (r e <C 1) 

The upper analysis obtains angular resolution of each EAS event individually. In case 
that there are many EAS events with different local coordinates which have contributions 
in the signature of each investigated source. Therefore at first angular errors of all of the 
accumulated EAS events in the galactic coordinates of the source were calculated, then 
the mean value of these angular errors was chosen as angular error of the source for 
the first step. Since all of the accumulated EAS events in the angular error region have 
contributions in the source signature, so the previous calculations were repeated for all of 
the accumulated EAS events in a circular region with the center of the source and radius 
of r e . Finally the mean value of these EAS angular errors is calculated as angular error 
of each source in galactic coordinates. Since the side distances of the array is different 
in El and E2, angular errors in these two experiments are different. So for calculation 
of the final result for each source these angular errors calculated separately for El and 
E2 and was weighted with the number of refined EAS events in the related experiment. 
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The final angular errors of investigated sources (r e ) are shown in Table 1. Since these 
angular error radii have a little fluctuations respect to their mean, therefore we sampled 
over 1 and b with a step of 5 degrees from the FOV of the array and calculated these radii 
to find the mean and standard deviation. Therefore the mean value and the standard 
deviation of the angular error of the experiment was obtained from angular error of more 
than 1000 random points. With These steps we obtained f e — 4.35° ± 0.82° as the mean 
angular error of the experiment. 

3.5. Drawing exposure map and simulation of the experiment 

Because of the various exposures of the sky in time, there is a non-uniform distribution of 
EAS events in galactic coordinates. The variation in time exposures due to the altitude 
difference of different sources, and the observation of separate individual galactic regions 
during the sub-intervals within the long duration of the experiment were simulated. We 
have 166665 EAS events in our experiments, so we used monte carlo method for this 
simulation and we simulated 166665 random events. These random numbers were chosen 
with considerations of Fig.|3 From this figure is seen that <j> distribution is not isotropic 
and the thickness effect of the atmosphere is very important, these effects were considered 
in choosing of these random numbers. So in the procedure : 

— Zenith angle (z) was taken from 1° to 60°. 

— Azimuth angle (</>) was chosen from 1° to 360°. 

— Related random numbers of time were chosen with consideration of EAS event rate 
of the experiment. Meanwhile we considered duration of the experiment which was 
taken from start and stop times of each sub-experiment. 

With this procedure we obtained 2500 simulated map and obtained the map with mean 
number of simulated events per (1° x 1°) pixel with the accuracy of 0.001. Fig.0]shows the 
exposure map of the experiment. The event map in Fig. |3 reflects the uneven exposure 
of the experiment. 

3.6. Investigation of EGRET gam ma- ray sources and measurement of their 
statistical significance 

The energy range of the logged EAS events by the array is in the range of 40 to 10,000 
TeV. In this energy range distribution of cosmic-rays is completely isotropic and ho- 
mogeneous in the galaxy. After correcting for the exposure effects, we looked for ex- 
cess emission that could be from gamma-ray sources. We used third EGRET catalogue 
IjHartman et al. 1999(1 as a reference. But some of EGRET sources have not acceptable 
events in the FOV of our array. So we counted number of events, number of pixels and then 
calculated count per pixel related to each of these sources. We note that the mean count 
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per pixel in data map is 4.798. Of 151 EGRET sources only 123 of them have count per 
pixel of more than square root of 4.798 with 98 more than 1.5 times the square root of the 
mean. So we started our investigations on these 98 sources. A method of excess similar to 
the analysis adopted by the Tibet EAS array has been adopted HAmenomori et al~ 2002 
fe l2(J(JUf) . In the first step we divided the data map (Fig.|3J) to the exposure map (Fig. 
pixel by pixel, in the obtained map, approximately most of non zero pixels are around 1 
except probable source pixels and pixels with more fluctuations in the data map, which 
probably due to the smallness of the data set. For eliminating the fluctuated pixels we 
multiplied the new map to 4.798 as raw exposure corrected map. In this step we added 
counts of all pixels of the raw corrected map. The number must be very near to 166,665 
so with this restriction we obtained a lower limit 0.0750 for eliminating pixels with less 
count in the exposure map, and the final exposure corrected map was obtained which is 
shown in Fig. [5] 

The obtained map was fairly uniform in the FOV of our array in the galactic coordi- 
nates. Next we investigated the remaining faint inhomogeneity in the corrected map that 
could be conditionally attributed to the existence of gamma-ray sources. To estimate the 
significance of an individual source we obtained all corrected EAS events, N on , within a 
radius v^fefrom the source position. The number of pixels, n s , within this region was 
also found. The total number of background counts, N a ff, was found from the pixels 
that fall within an outer radius of 2r e and an inner radius V2r e from the source position. 
The number of background pixels, rib, was also calculated. The statistical significance of 
the source was obtained using the Li & Ma relation IjLi fc Ma 1 983). 

S= N ™~ aN °ff- , a =^ (13) 
y/N on + a 2 N off n b 

Distribution of statistical significance of these 98 sources is fitted on a gaussian function 
as follow : 

f(a)=a a e X p(- i<T 2 ^ )2 ) (14) 

which is shown in Fig. [B] Then for the investigation of the statistical significance dis- 
tribution, we chose 98,000 virtual random sources with similar conditions of the 98 
EGRET sources. Fig. shows a normal distribution with mean 0.044 and standard 
deviation 1.001. Procedure of the significance calculation of these virtual sources is as 
like as the 98 EGRET sources except a little difference. Basically the virtual sources have 
no signals so we used another relation for the calculation of their statistical significance 
IjLi & Ma 19830 . 

S= ^»-^// _ , (15) 
y/a(N on + N off ) ri h 
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3.7. Investigation of a probable displacement of the source signatures 

Our exposure corrected map has not any bright source signatures, so we used third 
EGRET catalogue as a reference for searching some sources in our energy range. But 
EGRET energy range is from 100 MeV to 30 GeV and our energy range is from 40 
TcV to 10,000 TeV. So to search for some sources in our data we ought to hope that 
these sources have had a spread spectrum at least from EGRET energy range to our 
energy range, like blazars, BL Lac objects, Flat-spectrum radio quazars or etc. Since 
usually these sources in different ranges of energies have not the same places exactly 
and they have a little displacement, we searched around these sources with one degree 
displacement. This displaced 1 and b are shown in Table 1 for each source. It means that 
around each source with statistical significance more than 1 we tried 8 (1° x l°)pixels 
around it and chose the location with highest statistical significance. 

4. Results 

4.1. Explanation of the Field of View (FOV) in galactic coordinates 

The rotation axis of the Earth passes near the star Polaris, angular difference between 
Polaris and the rotation axis is approximately 5 times smaller than our mean accuracy 
(f~ e ) in galactic coordinates. So in this analysis Polaris is considered as being on the 
rotation axis of the Earth. The longitude and latitude of Polaris in galactic coordinates 
are 123°, 17 and 27°, 28 respectively. The geographical latitude of Tehran is 35° 2V, 
so the angle between the zenith of the array and Polaris in Tehran is 55°, we selected 
events with zenith angles less than 60° for the analysis which is deduced from Fig. H^b) 
and therefore, Polaris and regions around it are observable only with High zenith 
EAS events. From Fig. it can be seen that the best observable region is from 

10° to 40° of zenith angles. In Fig. |21 is shown that Galactic longitudes smaller than 
I w 123° - (2 x 60° - (60 - 55)°) « 8° and larger than I m 123° + (2 x 60° - 5°) « 238° 
are less observable. In other words, given the location of the array these are two 
different observable regions in galactic coordinates. Galactic latitudes smaller than 
b w 27° - (60° - 5°) w -38° and larger than b » 27° + (60° - 5°) » 82° are less 
observable regions too. 

4.2. Comparison of observed sources of El and E2 

With the procedure which is mentioned in subsection 3.7. we searched for sources with 
statistical significance more than 1.5, and we found thirteen sources which five of them 
are more than 2. For avoiding from probable fluctuations we did another try. we searched 
these displaced sources in El and E2 separately. But in this stage because of the smallncss 
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of these data sets, specially in El we separated significance more than 1. Therefore 10 
sources remained which have statistical significance more than 1 in El and E2 and more 
than 1.5 in the sum, which are shown in Table 1. Fortunately five of these sources are 
AGNs, one is probably AGN and four of them are unidentified sources. This is in case 
that from 271 source of 3rd EGRET catalogue sources only 66 are AGNs. 

4.3. Distribution of observed shower events around most significant EGRET sources 

It seems that radial distribution of number of counts per pixel for each source naturally 
must be near to a gaussian distribution as a source signature, over a flat back ground. We 
separated eight regions with approximately the same number of pixels for each source. 
The first region is a circle with radius yjl/2r e . The second region is a ring with inner 
radius \flj2r e and outer radius ^j2/2r e and with this order we separated eight regions. 
Distribution of mean counts per pixel around 98,000 virtual random sources and 10 most 
significant EGRET sources is shown in Fig. These distributions fitted on a gaussian 
function over a flat distribution as follow: 

f(r e ) = a r + b r exp{-r 2 e /2c 2 r ). (16) 
5. Discussion and concluding remarks 

In Table 1 is seen that most significant excesses observed, are in the region 15° < z < 35°. 
This result is reasonable because these angles are in favorable locations in the sky and 
have considerably more data from this region. In addition our data have a few counts in 
some parts of Fig. [3] and we were mandated to eliminate some source candidates from 
our list. To increase the statistical significance of our results and investigation of more 
sources, we have to accumulate more data to have a map with less fluctuations. 
There has been a considerable effort worldwide to detect gamma-ray sources via the 
EAS technique. From a variety of arguments we suspect that some, if not many, of the 
EGRET sources would be detectable at very high energies. In this work, we are lim- 
ited to a discussion of a few sources with relatively small statistical significance. Our 
statistical significance are not in a detection limit with confident, we studied this pro- 
cedure to guess some candidates in unidentified EGRET sources more than TeV range. 
For these sources listed in Table 1, we suspect that nine of them may be extra-galactic 
> 20°) QGehrels et al. 2000|) and only one is in galactic region (|6| < 20°) and this 
one is an AGN in the third EGRET catalogue list too. Four of our ten sources were 
investigated before with CASA-MIA (JOatanese et al. 1996J) and two of them were GEV 
EGRET sources l|Lamb & Macomb 1997jl . Therefore we might expect that as many as 
four of these unidentified sources could indeed be emitters at high energy and might be 
AGNs. 

Some of our observed sources overlap one another Fig. EI so a complete and accurate anal- 
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ysis procedure would incorporate the maximum likelihood method (Mat tox et al. 1996|) . 
We must also emphasize that our experiment can not distinguish between gamma-ray 
and cosmic-ray initiated air showers, and so we used the excess method to carry out a 
search for very high energy gamma-ray emission. After the analysis we understood that 
the record times per second of our computer is very important and we have to increase 
this record rate to decrease angular error radius of observable sources. In our future site 



at 2600 m a.s.l. (see http://sina.sharif.edu/~observatory/ 1, we are constructing under 



ground tunnels which will provide us with ample space to deploy muon detectors. The 
detection of muons in air showers should be provide a powerful away to discriminate 
between cosmic-ray and gamma-ray air showers. 
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Table 1. Observed EGRET 3rd catalogue sources by our array. Id and hd are dis- 
placed galactic coordinates, <?ei, &ei and Otot are the statistical significance related 
to the first experiment, the second experiment and sum of both, r e is error angular 
radius, z is mean amount of zenith angles of EASs related to each source and 'Flux' 
is number of EAS events related to each source, tl is AGNs which is investigated by 



CAS A- MI A before l|Catanese et al. 1996jl . t2 is sources with energy more than 1 GeV 
l|Lamb & Macomb 1997| . Meanwhile sources number 5 and 7 are 'Mrk 421' and '4C 
+10.45' respectively 
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Fig. 1. Experimental set up and electronic circuits. 
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Fig. 2. Local coordinates distributions of, (a) azimuth '</>' and (b) zenith V angles of 
logged EAS events. 
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Fig. 4. Exposure map of simulated events in l°x 1° bins base on the general parameters 
of total distribution of EAS events in galactic coordinates. 
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Fig. 5. Corrected exposure map which is extracted from division of the data map (Fig.[3J) 
to the exposure map (Fig. 0J pixel by pixel. 
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Fig. 6. Distribution of frequency of (a) 98,000 virtual random sources and (b) 98 EGRET 
sources in the FOV of our array in galactic coordinates versus their statistical significance. 
a<7, b a and c CT are described in eq. (14) 
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Fig. 7. Distribution of mean count per pixel of (a) 98,000 virtual random sources and (b) 
10 EGRET sources of Table 1 versus error radial distance from the centre of the related 
sources. 
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Fig. 8. Map of EGRET sources with statistical significance more than 1.5 in galactic 
coordinates. The numbered sources were described in Table 1. 
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